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Abstract: In this work we study the contribution of magnetic fields to the Sunyaev 
Zeldovich (SZ) effect in the intracluster medium. In particular we calculate the SZ angular 
power spectrum and the central temperature decrement. The effect of magnetic fields is 
included in the hydrostatic equilibrium equation by splitting the Lorentz force into two 
terms - one being the force due to magnetic pressure which acts outwards and the other 
being magnetic tension which acts inwards. A perturbative approach is adopted to solve 
for the gas density profile for weak magnetic fields (< 4/xG). This leads to an enhancement 
of the gas density in the central regions for nearly radial magnetic field configurations. 
Previous works had considered the force due to magnetic pressure alone which is the case 
only for a special set of field configurations. However, we see that there exists possible 
sets of configurations of ICM magnetic fields where the force due to magnetic tension will 
dominate. Subsequently, this effect is extrapolated for typical field strengths (~ 10/iG) and 
scaling arguments are used to estimate the angular power due to secondary anisotropics at 
cluster scales. In particular we find that it is possible to explain the excess power reported 
by CMB experiments like CBI, BIMA, ACBAR at i > 2000 with ag ~ 0.8 (WMAP 5 year 
data) for typical cluster magnetic fields. In addition we also see that the magnetic field 
effect on the SZ temperature decrement is more pronounced for low mass clusters ((T) ~ 
2 keV). Future SZ detections of low mass clusters at few arc second resolution will be able 
to probe this effect more precisely. Thus, it will be instructive to explore the implications 
of this model in greater detail in future works. 
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1. Introduction 

It has been known that the intracluster medium has magnetic fields of micro-gauss strength. 
They affect the evolution of galaxies [5], contribute significantly to the total pressure of 
interstellar gas, are essential for the onset of star formation [59], and control diffusion, 
confinement and evolution of cosmic rays in the intracluster medium (ICM) [32] . In clusters 
of galaxies, magnetic fields may play also a critical role in regulating heat conduction (e.g., 
[11, 42]), and may also govern and trace cluster formation and evolution. 

Magnetic fields in the intra-cluster medium have been inferred in various manners 
using diagnostics such as radio synchrotron relics within clusters, inverse Compton X-ray 
emissions from clusters, Faraday rotation measures of polarized radio sources within or 
behind chistcrs and cluster cold fronts in X-ray ([14, 8]). From these observations it can be 
inferred that the medium within most clusters is magnetized with typical field strengths 
at a few iiG level distributed throughout the cluster scale. In the cores of "cooling flow 
clusters" ([21, 64]) and also in cold fronts [65], the field strengths may reach 10//G-40//G 
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and could be dynamically important. It is thus essential to study the effect of an intra- 
cluster magnetic field on the gas density distribution in the ICM. 

A very important observational probe of the cluster gas density distribution is the 
thermal Sunyaev Zeldovich effect [62, 63]. This effect occurs because of the re-scattering of 
primary CMB photons with the hot electrons in the ICM resulting in secondary anisotropics 
in the CMB at small angular scales (for a detailed review see [6]). At 30 GHz, anisotropics 
from the thermal SZ effect are expected to dominate over the primary CMB fluctuations 
for multipoles i < 2500. The angular power of SZ fluctuation depends sensitively on the 
integrated cluster abundance and cluster gas distribution (for details see [2]). The SZ effect 
in clusters depends directly on the density and temperature profile of the ICM. For any 
particular cluster, it is quantified through the y-parameter which is essentially the temper- 
ature decrement along a line of sight through the cluster. This is therefore an important 
observational probe and can be used to constrain/detect any additional parameter that 
affects the density of the ICM. In the present work, we study the effect of intracluster mag- 
netic fields on the gas density profile of the ICM and compute the central SZ decrement for 
different cluster masses as well as the angular power spectrum . To study magnetic field 
effects in the ICM, earlier authors have incorporated it by introducing the magnetic pres- 
sure in the hydrostatic equilibrium condition ([28, 15, 67]). Although this is an important 
effect, there is, in addition to the pressure, a contribution arising from the tension force 
due to the field. We show in this paper, that the most general treatment involves using 
both terms arising out of the Lorentz force due to the field. In general, the detailed form of 
the magnetic field configuration is unknown and moreover it can vary from cluster to clus- 
ter. It is only in the specific case of isotropic configurations, that the magnetic pressure is 
the only contribution. However, there could be other plausible configurations of magnetic 
fields in clusters where both the pressure and the tension force contribute. In this paper, 
we classify the different field configurations and study the problem for a near-radial field . 

The paper is organized as follows: Section 2 presents the formulation of the problem by 
setting up the hydrostatic equilibrium equation incorporating the magnetic field. Sections 
3 and 4 present the mathematical formulation of the y-parameter and the CMB angular 
power spectrum respectively. Section 5 presents the results for the density profiles and 
SZ observables which are got by solving the hydrostatic equilibrium equation and finally 
section 6 summarizes the results and concluding remarks. 

2. Formulation 

In this section, we present the method of incorporating the effect of a generic tangled 
magnetic field for studying its effect on the ICM gas. We assume that the ICM is in 
hydrostatic equilibrium in the presence of magnetic field whose effect is incorporated by 
using the Lorentz force. The typically smooth morphology of the X-ray emission from 
the hot intra-cluster medium leads naturally to the hypothesis that the gas is in near 
equilibrium, stratified along isopotential surfaces in a mildly evolving distribution of dark 
matter, gas and galaxies. This suggests that the assumption of hydrostatic equilibrium for 
such relaxed clusters is mostly justified. The dark matter mass profile and the temperature 
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profile are specified by their 'universal' forms as described in the following sections. The 
radial profile of the magnetic field is also implicitly specified by its dependence on the 
density profile. Given these profiles, the gas density is finally solved perturbatively around 
the default profile (i.e the profile in the absence of magnetic field). The perturbation 
method is only applicable for small magnetic field values (in this case, it can be used for 
values of the central field upto ~ 4:fiG) . 

2.1 Dark matter density profile 

The gravitational potential within clusters is mainly determined by dark matter. The dark 
matter density profile, Pdmir) suggested by many high resolution N - body simulations is 
well described by the NFW profile ([43, 44]) 

Pdm = A2/dm(a;) = (2.1) 

Here, x = r/r^ is the radius in unit of core radius and is a normalization factor which 
represents the characteristic density at , r = r^. Since the dark matter density profile 
is self-similar, the dark matter mass profile is also self-similar. So, the dark matter mass 
enclosed within a radius r is 

M{< r) = 47rp,rfm(r/r,) (2.2) 
where, m{x) is a non-dimensional mass profile given by 



X 

m{x) = / duu^ya^iu) = ln(l + x) - ; (2.3) 

Jo U + ^1 

The definition of the virial radius, i?vir, is the radius within which the total dark matter 
mass is confined, i.e., M^ir = M(< c), where 

c=— (2.4) 

r 

' s 

is a dimensionless parameter called the 'concentration parameter'. Evaluating equation 
(2.2) at the virial radius, the normalization factor, p^, is fixed at; 



The virial radius, ii^ir {M^ir, z) is calculated with the spherical collapse model [45] 



^ir = 



- 1/3 




1/3 




_47rpsm(c) 





(2.6) 



where the second equality comes from evaluating R^i^ from equation (2.5). Here A^{z) is 
the spherical over-density of the virialized halo within R^^^ at z, in units of the critical 
density of the universe at z, Pc{z). Following [29], we assume a value A^{z = 0) = 100 for 
a cosmological model with = 0.29 and = 0.71. 
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We follow [7] in adopting the approximation for c as a function of the virial mass of 
the cluster. They give the median values of 'c' and also the la deviations: 

/ M- \-0.13 

- = ^( i,5x10.3;-imJ (") 

with K = 9 reproducing the best-fit and K = 13.5 and K = 5.8 reproducing the +la 
and the —la values in the concentration parameter. These values of the concentration 
parameter are also consistent with the findings of [60] . 

The above set of equations specify the dark matter density profile of a particular mass 
cluster. 

2.2 Universal temperature profile 

The "universal temperature profile" used for our calculation [34] is (normalized by the 
emission- weighted temperature): 

I) = (iT^ ^'-'^ 

where (T) is the emission-weighted temperature of the cluster, Tq = 1.33, Ux = R^i^/l.b, 
and 6 = 1.6 on the radial range (0.04-1.0) i^^ir- To determine the emission-weighted 

temperature from the cluster mass, we use a relation that arises from adiabatic evolution 
of the gas in cluster. [1] have shown that the observations of [23] of Mr^^f„-{T) relation 
in clusters can be understood from gravitational processes alone. We therefore use this 
empirical relation (M5oo-(r)) derived by [24]: 



M500 = (2.64iO;l)10i3 ( ^-^^^-^ ^^.9) 

where fcb is the Boltzmann constant and M500 has been calculated self-consistently by taking 
the total mass within the radius where the over-density is (5 > 500. 

2.3 Hydrostatic equilibrium of gas in presence of magnetic field 

In this section, we determine the density profile of the gas in the ICM given the universal 
temperature profile and assuming that the gas is in hydrostatic equilibrium (HE) with the 
background dark matter in the presence of a tangled magnetic field specified by a profile 
B{r). 

2.3.1 Density profile in the absence of magnetic field 

The HE equation for the gas in the absence of a magnetic field has the well-known form: 

1 j^^_GM(<0 

Pg{r) dr 

where G is the gravitational constant, /9g(r) and -Pg(r) are the density and pressure profiles 
respectively. The pressure is related to the density through the equation of state: 

^ p{r)kBT{r) 
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Here, M(< r) is the mass enclosed within radius r and /x and nip denote the mean molecular 
weight (// = 0.59) and the proton mass. T{r) is the universal temperature profile of the 
gas. M{< r) is mainly determined by the dark matter mass profile. The solution p^{r) to 
the above HE equation is referred to as the default density profile. Using the equation of 
state, the HE equation can be recast as: 

kB fdlnpg ^ dlnT{r) \ ^ GM{< r) ^2.12) 



pnip \ dr dr J r^T(r) 

The solution Pg(r) to the above can then be expressed as : 



Ps(0) T(r) ' V kB J„ r^T(r) 

The central density Pg (0) is fixed by the constraint that the gas mass within the cluster is a 
'universal' fraction of the total mass of the cluster. Thus, the gas fraction /g^s = M^^/M^^^^i 
within the virial radius is taken to be universal and equal to 0.105, as recently found by 
Ettori (2003) for a sample of low and high redshift clusters. The gas mass can be expressed 
as: 

Mg,, = 47r / drrVg(r-) (2.14) 
Jo 

Since the total gas mass is negligible compared to the dark matter, M^^^^x ^ M^^, mass in 
dark matter, and therefore /ga^ ~ M^^/M^^. 

2.3.2 Density profile in the presence of magnetic field 

The effect of magnetic fields is now taken into account by introducing an additional term 
corresponding to the radial component of the Lorentz force due to the field. The Lorentz 
force due to a tangled magnetic field B{f) is given by 

F{f) = B{f) X (V X B{f)) (2.15) 

The Lorentz force depends on the magnetic field configuration which is not known and 
moreover its detailed form could vary from cluster to cluster. However since we are looking 
at only the radial profiles, we incorporate the effect of magnetic field by considering the 
angle-averaged radial Lorentz force {F^). We then simplify the angle-averaged Lorentz force 
by explicitly writing it down in terms of the angular averages of the field components (i.e in 
terms of {Bf), {Bg), {B^)). This method has been used in the context of studying magnetic 
field effects on the spherical collapse of non-rotating low-mass gas clouds [12] and also in 
the context of spherical accretion onto a black- hole in the presence of magnetic field [58] . 

It is important to note here, that the Lorentz force due to magnetic field consists of 
two terms, one term which is a pure gradient V-B^, known as magnetic pressure and the 
other one being {B.V)B known as magnetic tension. Although these two terms contain 
the gradient operator V, they are not identical. The magnetic pressure is a pure gradient 
dB"^ /dr^ and hence is directly tied to the gradient of the underlying gas density distribution. 
So for instance, if the gas density profile is nearly flat, then magnetic pressure has a 
negligible contribution. On the other hand, the magnetic tension is not a pure gradient. 
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The magnetic tension component along the radial direction can be written as: (B.'V)Br = 
V.(-B-Br)) (using the fact that V.-B = 0). Hence this tension component is a divergence and 
hence consists of two sub-terms, one of which is like dB'^/dr (pressure-like term) and the 
other is like B'^/r. This non-pressure like term also known as "hoop stress" acts opposite 
to the pressure gradient and increases towards the center. Refer to Eq All and A12. We 
use eq A12 for the derivation of dcnsitj^ profiles. We have also included the very special 
single case of equal field components also for reference. It is only in this case that the 
Lorentz force has a pressure-like term. But in our work we consider generic cases in which 
the hoop stress contributes. 

We classify the different sets of field configurations in terms of the relative strengths 
of the radial and transverse correlations (For details refer to the Appendix). Previous 
studies have mainly focused on the isotropic configuration in which the relative strengths 
of the three correlations are equal. In this case it can be shown, that the angle- averaged 
Lorentz force acts like a magnetic pressure. However this is only a special configuration. 
We consider a generic set of configurations in which the radial component of the magnetic 
field is auto-correlated , while the correlation between the transverse(polar) components is 
negligibly small. 

To begin with, the HE equation in the presence of tangled magnetic fields has the 
form: 

1 dPs{r) GM{<r) 



Pb (r) dr 

In this case too, the equation of state is specified by 



+ (F.) (2.16) 



where Pb{t), P-nif) TB{r) are the modified density, pressure and temperature profiles 
respectively in the presence of the magnetic field. 

In order to attempt to solve the equation for the modified density profile Pb(^) we 
make some simplifications. Firstly we assume that the field-strength is not very large so 
that a perturbative solution around the default density profile p,{r) can be attempted. 
For this to be true, the energy density in the magnetic field should be small compared 
to the gravitational energy density corresponding to the default gas density. The next 
thing to note is that the change in the default temperature profile will be negligible (i.e 
T-Bir) = Tg{r)) since we do not consider any additional energy transfer to the gas from 
the magnetic field. We use the universal temperature profile as an input to deduce the 
density profile since its a good fit to observations. Since we use a perturbative scheme, the 
change in the temperature profile due to magnetic field contributes at only higher order in 
perturbation theory and hence can be neglected. Under these assumptions we can combine 
the two equations (2.10) and (2.16) to get: 
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In this case, the gas fraction f^^^ = M^^^/M^^^^^i is considered to be 'universal' too and equal 
to 0.105. The upper limit a in the above integral is determined from the relation: 



In order to compute the modified profile psir), it remains to specify the radial profile 
of the angle-averaged Lorentz force {F^{r)). As shown in the Appendix, {F^), in turn, can 
be expressed in terms of the angle-averages of the squared magnetic field components viz. 
{B"^), (Bg) and (-B^). Thus the form for F^{r) depends on the relation between the different 
component profiles. 

In the present work, we compute the results for field configurations which are nearly 
radial, i.e B^ = B^. For such configurations, the form for F^{r) is given in eq(A12). It 
thus remains to specify only the field strength profile B(r) and this is done below. The 
magnetic field profile is assumed to be of a parametric form, 



Here, Pg{z = 0) is the average cosmological density of the gas at z = 0, B^ is measured 
in and a is a parameter > 2/3. The above form is suggested by the rotation measure 
(RM)-X-ray correlation inferred in simulations/observations, (see, e.g., [8, 20]). In this 
paper, we assume the value a = 0.9 for our calculations. ([20]). 

3. Compton y-parameter and central SZ decrement 

The temperature decrement of CMB due to the SZ effect is directly proportional to the 
Compton parameter (y). For a spherically symmetric cluster, the Compton parameter is 
given by 



where ctt is the Thomson cross-section, and Pe{r) = ne{r)k\^Te{r) is the electron pressure of 
the ICM, where ne(r) = 0.875(/9gas/mp) is the electron number density, fcb is the Boltzmann 
constant, and Tc(r) is the electron temperature. The integral is performed along the line-of 
-sight [l) through the cluster and the upper limit of the integral {+R) is the extent of the 
cluster along any particular line-of-sight. We do not include the effects of beam size in 
calculating the y parameter. This approximation is justified by the fact that the pressure 
profiles are relatively fiat in the inner region. The variation of pressure integrated along 
the line-of-sight as a function of the projected radius is even flatter thus providing more 
justification for the above approximation. 

The angular temperature profile projected on the sky due to SZ effect, Ar(^)/TcMB 
is given in terms of the Compton parameter in equation (3.1) 




(2.19) 




(2.20) 




(3.1) 



Ar(e) 



= 9{x)y{e) 



(3.2) 



-7- 



where g{x) = xcoth(.T/2)-4, x = hi^ /k^TcMB, ^cmb = 2.728 [25]. In the Rayleigh-Jeans 
approximation, g(x) ~ —2. We only evaluate "central" SZ decrement from the pressure 
profiles of our models. In this case, the integral in equation (3.1) reduces to 

yo = / Pe{r)dr (3.3) 

ii'-eC Jo 

In the Rayleigh-Jeans part of the CMB spectrum, the deviation from the black-body spec- 
trum results in a decrement of the CMB temperature, 

Ar^wO«-5.5yiC (3.4) 

We use the pressure profiles resulting from our model to calculate the central SZ decrement 
in the temperature of the CMB. 



4. Angular power spectrum 

The angular two-point correlation function of the SZ temperature distribution in the sky 
is conventionally expanded into the Legcndrc polynomials: 

/ AT AT 



(n)-^(n + 0)) = i- J^(2^ + l)C,P,(cos^) (4.1) 



\2CMB TCMB / 47r 

(P) (C) (P) 
Since we consider discrete sources, we can write = Q + w \ where Q is the 

contribution from the Poisson noise and c}*^^ is the correlation among clusters ([45], § 41). 

We define the frequency independent part in the power spectrum as C^^^^ = Ci/g^{x). 

The integral expression of C^*^^-* can be derived following [16] as 

where V{z) is the co-moving volume and yt is the angular Fourier transform of y{9) given 

by 

= 27r y" y{9)U{t + l/2)e]ede, (4.3) 

where Jq is the Bessel function of the first kind of the integral order 0. In equation (4.2), 
Zdec is the redshift of photon decoupling and dn/dM is the mass function of clusters which 
is computed in the Press-Schechter formalism [51]. The mass function has been computed 
using the power spectrum for a A-CDM model with normalization of (jg = 0.8. We choose 
Mnin = 5 X IO^^Mq and M^^ = 2 x lO^^M© and integrate till redshift of z = 5 instead 
of Zdoc- This is done because the integral in equation (4.2) is found to be insensitive to 
the upper limit in redshift beyond z = 4, the reason being that the mass function is 
exponentially suppressed beyond that redshift in this mass range 
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5. Results Discussion 



In this section, we investigate the effects of the magnetice field with nearly radial config- 
uration on the density profiles, the central SZ decrement and on the SZ angular power 
spectrum. We also present a discussion on the other possible confiurations of the cluster 
magnetic field and the Lorentz force arising out of thes configurations. 




Figure 1: Pressure (Pmag) and tension 

(Tmag) forccs arising out of the mag- 
netic Lorentz force as a function of the 
scaled distance for different magnetic 
field configurations like transverse (tr), 
isotropic (iso), anisotropic (aniso) and 
nearly radial (rad). The pressure and 
tension terms are computed with arbi- 
trary normalization 



4 




Figure 2: Angle averaged Lorentz force 

calculated for different magnetic field 
configurations. All the nunerical factors 
have been normalized to unity 



5.1 Magnetic field configurations and Lorentz force 

The Lorentz force can be split into two terms : the magnetic pressure (a d{B'^)/dr) 
which depends on the gradient of the field strength and the magnetic tension (oc B^/r) 
which depends on the field strength (as shown in detail in the Appendix). In figure (1), 
we plot the radial dependence of the magnetic pressure and the magnetic tension for dif- 
ferent field configurations viz. radial {{Bt^) = 0), transverse ((-Br^) = 0), anisotropic 
{{Bt"^) = Q.b{B^))axid isotropic ((5^^) = {Be^) = (B^) = 0.33(52) y ^^^^^ (^2), we 
plot the resultant Lorentz force (sum of pressure and tension) for the same set of config- 
urations. To demonstrate the behaviour of magnetic pressure and tension forces as well 
as the angle-averaged radial Lorentz force as a function of scaled radius of the cluster, we 
have normalized all the constants to unity in both the figures. Wc also use the convention 
that positive values correspond to force directed radially outwards whereas negative values 
correspond to force directed radially inwards while plotting. As seen in figure (1), for the 
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radial configuration, the tension force dominates over tlie pressure force and the net force 
is directed radially inwards as we approach the central region (r < 0.2r2oo) of the cluster. 
This causes an enhancement in the gas density relative to the default state. This is different 
to the manner in which the effect of the magnetic Lorentz force is usually incorporated in 
the hydrostatic equilibrium equation for the cluster gas. 

For the transverse configuration, the tension force dominates over pressure force and 
is directed outwards as a result of which, the net force is directed outwards as we approach 
the central region. This will possibly cause a decrement in the gas density relative to the 
default case. In such models, the magnetic field can act as a heating source in the ICM. 

In the anisotropic case, the pressure force vanishes whereas the tension force is directed 
inwards and thus the net force is directed inwards. In the isotropic configuration, the 
tension force vanishes whereas the pressure force is directed inwards and thus the net force 
is directed inwards. In both these cases, the gas tends to fall towards the central region 
thus causing an enhancement in gas density relative to the default scenario. 

0.1 E . • j 




Scaled Radius r/rggo 

Figure 3: Density versus scaled radius r/raoo for cluster masses 5 x IO^^Mq and 1 x lO^^M© 
5.2 Density profiles 

Using the NFW profile for dark matter and the 'universal' temperature profile, the hydro- 
static equilibrium equation is solved to get the density profiles for different magnetic field 
strengths for a given cluster mass. We consider the range of cluster masses 5 x 10^^ — 1 x 

In figure (3), the density profile (in units of particlcs/cc) is plotted as a function of 
the scaled radius r/raoo for a clusters of mass 5 X lO^-'^M© and 1 x IO^^Mq for different 
values of field strengths specified by 5*. As seen in figure (3), the presence of magnetic 
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Figure 4: Fractional change in density due 
to magnetic fields of strength 3 fiG and 4/iG 
versus scaled radius r/r 200 for a cluster mass 
5 X IO^^Mq 



Figure 5: Fractional change in density due 
to magnetic fields of strength 3fiG and 4/iG 
versus scaled radius r/r^oo for a cluster mass 
1 X IQ^^Mq 



field causes an enhancement in the gas density compared to the default (i.e no magnetic 
field) value for a particular cluster. The enhancement is large in the central regions and 
decreases outwards. There is a crossover from enhancement to depletion which occurs at 
TcTo = 0.2r2oo- At distances greater than the crossover radius there is a depletion in the gas 
density. This occurs because of the fact that the Lorentz force contributes through two 
opposing forces viz. the magnetic tension and magnetic pressure. The magnetic tension 
dominates in the central regions whereas magnetic pressure dominates in the outer regions. 
The magnetic tension acts inward and hence pushes the gas inwards increasing the density 
whereas in the outer regions (where magnetic pressure dominates) the magnetic force acts 
outward and pushes the gas outwards, depleting the density. 

In figures 4 and 5, the fractional percentage density change is plotted for masses 
5 X IQ-'^'^Mq and 1 x IO^'^^Mq respectively. There is nearly a 70% increase in the density 
towards the central region of the cluster of mass 5 x lO^^M© for a magnetic field strength of 
4/xG. For a higher cluster mass this fractional change in density is lower as seen in Figure 3. 
The extent of density enhancement for a given magnetic field strength thus shows a trend 
with cluster mass, the enhancement decreasing with increasing mass. In other words, the 
effect of a magnetic field of given strength on the gas density is higher for low-mass clusters 
and groups. 

This enhancement of density in the central regions has important implications for the 
Sunyaev Zeldovich effect since the temperature decrement for an individual cluster depends 
directly on the gas density. This in turn also gets reflected in the CMB angular spectrum 
on small scales or high multipoles. We discuss these issues in the following sections. 



- 11 - 




1 2 3 4 5 6 7 8 

(T) keV 



Figure 6: Observed and predicted relations of central SZ decrement {AT^yjo (/xK)) versus emission- 
weighted temperature ((T) (keV)) in clusters. The solid line corresponds to magnetic field strength 
= 4/iG, the dotted line corresponds to = S/iG and the dot-dashed line corresponds to 
S* = 0/LtG. Data points have been taken from Zhang & Wu 2000, McCarthy et al. 2003 and 
references therein and Lieu, Mattiz, & Zhang 2006 

5.3 Central SZ decrement 

Using the density profiles, the pressure profiles are evaluated and the resulting central SZ 
decrement is computed as a function of the cluster mass. In figure (6), the central SZ 
decrement is plotted as a function of the emission-weighted temperatures of the cluster . 
For comparison, we have also plotted the data for the central SZ decrement for a sample 
of clusters taken from [68, 39] and references therein and observed with WMAP [33]. 

From figure (6), it can be seen that the temperature decrement is enhanced in the 
presence of magnetic fields and increases with the emission weighted temperature (and 
hence the cluster mass). The effect of the magnetic field is more pronounced for low mass 
clusters. In particular, magnetic fields of strength ~3/LtG cause nearly a 25 % increase in 
the SZ decrement for a cluster with emission weighted temperature ~2 keV corresponding 
to a virial mass of 1 X IO^^Mq. Thus from precise SZ observations of low mass clusters 
and galaxy groups it might be possible to detect or place constraints on the field strength 
in these structures. 

Currently proposed SZ experiments are expected to make wide angle surveys of the 
CMB sky with a resolving power of upto few arc minutes . Experiments such as South Pole 
Telescope (SPT) and Atacama Cosmology Telescope (ACT) can survey several hundreds of 

square degrees upto a mass limit 2x10"*^^ Mq with high sensitivity at arc minute resolution 
[56] . In particular, the SPT survey will cover 4000 deg^ of sky with lO/^K sensitivity per 
1 arc minute pixel at 150 GHz and is expected to yield nearly 10000 clusters with masses 
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greater than 2xlO"'^^M0 [9]. The ACT survey is expected to map 200 square degrees of 
CMB sky in 3 frequency bands of 145 GHz, 220 GHz and 265 GHz at arc minute resolution 
reaching sensitivity levels of 2//K per pixel [37]. However to detect the CMB temperature 
decrement in the central regions it is necessary to be able to resolve features at the level of a 
few arc seconds. This is because the cluster virial radius (~1 Mpc) corresponds to angular 
scale of arc minutes whereas the central regions (~ 100 kpc) will typically correspond to 
an angular scale of few arc seconds. There arc upcoming proposals which will probe such 
scales. In particular, the Cornell-Caltech Atacama Telescope (CCAT) [26] is proposed to do 
high angular resolution follow-up observations of clusters which have already been detected 
in some of the above mentioned wide-area surveys. The 150 GHz camera on CCAT would 
be sensitive at 310 /iK s^/^ with an angular resolution of nearly 3 arc seconds per pixel. 
These experiments will thus provide a statistically large sample of clusters using which it 
would be possible to detect or put proper constraints on the ICM magnetic field strengths. 



5.4 SZ power spectrum 
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Figure 7: Poisson contribution to the Figure 8: Fractional change in angu- 

angular power spectrum of the SZ flue- lar SZ power for different magnetic fleld 

tuations plotted as a function of f strengths plotted as a function of i 

The results for the SZ angular power spectrum are presented in figure (7). The power 
spectra are plotted for different magnetic field strengths and the power spectra for the 
default density and temperature profiles is also plotted for comparison. It can be seen that 
the power peaks close to I = 3300 for a field strength of 3)uG and the peak shifts to higher 
i values as the field strength increases. This means strong magnetic fields have a dominant 
effect on small scales. The typical value of the SZ power for a field strength of 3/iG is 
3.2fiK^. The SZ power increases with increasing magnetic field strength. In the figure, we 
have plotted the results for magnetic field strengths upto 6fiG. The perturbative estimate 
is applicable only for field strengths upto 3/LiG. The SZ power increases with increasing 
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magnetic field strength as can be seen in the figure. We can extrapolate this trend to 
estimate the CMB power at larger values of field strengths (~ the order of few-tens of 
micro-Gauss, as has been observed in clusters) by a simple qualitative scaling argument. 
Since the magnetic field depends on density as 5 oc p^/^, and the y-parameter (which is the 
fractional temperature change) is directly proportional to density AT/T oc p, we deduce 
that the CMB angular power which scales as temperature squared will scale with magnetic 
field as Ci oc i?^. From this we can see that for magnetic fields of strength ~ 15-20;uG, 
the CMB power will increase by a factor of nearly 125-300 times compared to the power 
estimated for 3 /xG. This means that the linearly extrapolated estimate of CMB power 
for a 15/xG field will be ~400/xiC^. In figure (7), C( for field strengths 5/xG and 6/LtG are 
plotted using this scaling estimate. 

In figure (7), we have plotted the percentage change in the angular power of CMB 
compared to the default case. It can be seen from the figure that a 3/xG field induces a 
20% increase in the CMB power whereas a 4/iG induces a 40 % change. 

Observations by CBI [52] and BIMA [18] report an excess in CMB anisotropy at 
£ > 2000. The primary anisotropies at these scales are damped out and the cause of 
this excess is attributed to the SZ power spectrum from galaxy clusters. CBI 5 year data 
collected at 30 GHz [61] report the level of anisotropics at £ ~ 3000 to be 380 fiK^ with a 
Gaussian error of 117 fiK^. An excess at nearly 1 a level is also reported by the ACBAR 
experiment [54] at 150 GHz at these scales. Attributing this excess to the SZ effect alone 
(without the inclusion of magnetic fields) would require the normalization of matter power 
spectrum (jg ~ 1. WMAP 5 year results [30] however have ruled out this value and fixed 
as ~ 0.8. As a result, this excess cannot be explained by the ordinary SZ effect alone and 
additional possible effects which might contribute have to be investigated (see for e.g [35]). 
However as we discussed above, if magnetic fields are present generically at the level of 
tens of fiG, then it might be possible to reconcile with this cTg ~ 0.8. 

Recent observations conducted by the Sunyaev Zeldovich Array (Sharp et al. 2009) do 
not report any excess and estimate the constraint on the level of secondary anisotropies 
to be li'^Q^fxK'^ . APEX-SZ has measured the CMB angular power at a frequency of 150 
GHz in the range of 3000 < £ < 10000 and report anisotropies at the level 33^24^-^^ 
at effective £ ^ 5000 [53]. These experiments detect much less power compared to CBI, 
BIMA etc. The CMB anisotropy power at this value would then limit the generic magnetic 
field strength to ~ 7 — 8/LtG. Several upcoming CMB experiments, such as the Atacama 
Cosmology Telescope(ACT), Planck etc claim to reach upto the cosmic variance limit and 
measure the SZ power to an accuracy of nearly 1 % (see the science white paper [41]). 

The central SZ decrement and the angular power spectrum depend strongly on the 
magnetic field configuration as well. In this work we have explored the effects of a nearly 
radial magnetic field. However, there could be other possibilities in which the effects could 
be different and might even be opposite to our calculations. The extreme case where the 
tension force goes to zero and we are left with only the pressure term is the isotropic 
configuration which has been addressed in earlier works by several authors [15, 28]. Here 
there is a suppression in the central SZ decrement and thus in the angular power, since the 
magnetic field will act as a heating source and push the gas out of the central regions. In the 
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other two cases, the effect of the magnetic field on central SZ decrement and angular power 
will vary according to the geometry of the field (as explained in the the figures (1) and 2. 
Thus this can be possibly used as an indirect probe of the magnetic field configurations in 
the intracluster medium. 

Thus it is extremely important to model the effects of intracluster magnetic fields in 
addition to the standard effects in order to calculate the SZ power spectrum much more 
precisely for different field configurations. 

6. Conclusion 

In this work, we have studied the effect of an intracluster tangled magnetic field on the gas 
distribution in the ICM. In addition, we have also investigated the effect of cluster magnetic 
fields on the central SZ decrement for a range of cluster masses. We have also computed the 
CMBR angular spectrum for different field strengths . In contrast to previous studies, we 
incorporated the complete radial Lorcntz force due to the magnetic field in the hydrostatic 
equilibrium equations, thereby introducing the effect of magnetic tension as well in addition 
to magnetic pressure for a generic field configuration. It is only in the special case of an 
isotropic configuration that magnetic pressure is the sole contribution. However, realistic 
cases would involve a range of configurations. In particular we presented the results for 
a nearly radial magnetic field configuration. The results would be more pronounced for a 
configuration of magnetic fields in which the tension force is dominant. 

It would be interesting to look at all the other important baryonic effects present in 
the cluster gas like AGN heating, cosmic-ray heating and the like in addition to this effect 
of the magnetic field to determine the structure of the intracluster gas in more detail and 
to disentangle various contributions. 

Our results can thus be summarized as follows: 

1. For a nearly radial magnetic field, the ICM gas density shows an enhancement in 
the regions close to the center whereas there is depletion in the outer regions, the 
crossover scale being dependent on the cluster mass. 

2. The gas density enhancement /depletion is large for low-mass clusters. 

3. The central SZ decrement is enhanced compared to the default (i.e no magnetic 
field) case if magnetic field effects are included because of the enhancement in the 
gas density 

4. The CMB angular power is also enhanced in the presence of magnetic fields and future 
precise observations on these small scales can be used to constrain the strength of 
such fields and their configurations. 

5. For an isotropic field configuration, the effects can be opposite to what we have 
concluded in this piece of work, as pointed out by earlier authors. 
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A. Angle-averaged radial Lorentz force 

The Lorentz force due to a magnetic field B is given by 

^_ ^x(Vx^) 
An 

This can also be written as : 



(A.l) 



AttF = ^V{B^) - {B.V)B (A.2) 

In spherical polar coordinates, the Lorentz force can be written in terms of its components 
as F = fF^ + dFg + ^F^ In particular, the radial component of the Lorentz force Fj. can be 
expressed as: 

B, dB, B,d Be d BodB, 

^^^r = — ^^rr ^{rB^} — {rBe) + ^ (A.3) 

rsm^ 0(f> r or r or r ou 

Adding B^{V.B) to the R.H.S this can be further rewritten as: 

4nF^ = V.iBj)-^-lliB^) (A.4) 

Here, = B^^ + Bg^, denotes the transverse part of the expression. We now evaluate 
the angular average of F^. The angle averaged Lorentz force is defined as: 

{F.) = j -F^ (A.5) 

Taking the angular average of Eq (A.4), we get: 

47r(F,) = (V.(i?.^)) - ^ - \^^{{B^)) (A.6) 

To get a simplified expression for the above angular average, the angle average (V.(BBJ) 
remains to be evaluated. This can be done in the following manner. Let 

/(r) = (V.(S5,)) = j ^V.{BB^) (A.7) 

Multiplying both sides by r^dr and integrating we get, 

J drr'fir) = J ^V.(SSO (A.8) 
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Using Stokes theorem, the right hand side of the volume integral above can be expressed 
as a surface integral and hence, 

J drr'^fir) = r'^ J ^f.{BB,) (A.9) 

Now, differentiating both sides w.r.t r we get: 

fir) = -Mr^{B^)) (A.IO) 

Using Eq (A. 5), Eq (A.9) and the fact that B"^ = B"^ + Bq + B'^, the expression for the 
angle-averaged radial force can be written as 

MF.) = ~ m - 2{B,')] + '^^'^-/^^-'^ (A.11) 

The above expression for the angle-averaged radial Lorcntz force can be simplified for some 
special magnetic field configurations which we enumerate below: 

1. {Be^) = (B/) = : 

In this case, (i?^) = (-Br^)- Hence 

In this case, the simplified expression is. 
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